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Abstract 

A micromechanical model based on the theory of elasticity has been developed to study the 
configuration of self-assembled impurity nanostructures in high temperature superconducting 
YBa2Cu307_(5 films. With the calculated equilibrium strain and elastic energy of the impurity 
doped film, a phase diagram of lattice mismatches vs. elastic constants of the dopant was obtained 
for the energetically-preferred orientation of impurity nanorods. The calculation of the nanorod 
orientation and the film lattice deformation has yielded an excellent agreement with experimental 
measurements. 
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I. Introduction 

Self-organization of nanostructures in epitaxial films may provide a unique approach 
to design and tailor physical properties of nano-composite films by controlling the mor- 
phology of the nanostructures. Impurity doped high temperature superconducting (HTS) 
YBa2Cu307_5 (YBCO) film is an excellent example among other oxide nano-composites. In 
pursuing electronic and electric applications of HTS materials, extensive efforts have been 
taken to improve the critical current density ( Jc) by incorporating impurity nanostructures 
as strong magnetic vortex pinning centers in YBCO films HQ. Self-organized impurity 
nanostructures, such as BaZrOa (BZO) nanorod (NR) array oriented along the YBCO c 
axis, were found to dramatically improve with up to a fivefold increase in applied mag- 
netic field up to 5 T jsl- The strong correlated pinning provided by the aligned NRs also 
reduces the Jc anisotropy in the applied field direction. BZO has, however, a relatively large 
lattice mismatch with YBCO, resulting in a reduced critical temperature due to a seriously 
strained YBCO lattice upon BZO doping ^. While a strained lattice is essential to the 
formation of the nanostructures, a quantitative control of the strain at a microscopic scale 
could optimize physical properties of nano-composite films. It is thus desirable to have quan- 
titative criteria in selecting dopants for a given film matrix based on an understanding of the 
strain field and its role in the formation of the nanostructures. Such criteria are, however, 
not currently available. Moreover, the perfectly aligned NRs along the c axis in YBCO films 
is not ideal since the in-field Jc is not much improved or even reduced when the applied field 
is not along the c axis. A three-dimensional pinning landscape in the film could be a solution 
to this problem and some attempts have succeeded recently in generating splayed BZO NRs 
using either vicinal substrate jsl or secondary impurity nanoparticles j?]. In both cases, an 
additional lattice mismatch is introduced into the film matrix to change the strain field in 
the film. Theoretical understanding of the interplay of strains due to multiple mismatched 
lattices becomes essential to the control of the microstructure and physical properties of 
nano-composite films. Moreover, a better modeling of the strain field in the doped YBCO 
films could also help to understand the change of the superconducting properties of the films 
due to the doping [8|. 

Many theoretical studies have been made on self-organization of nanostructures in crystals 
jo], lot. The analysis based on the elastic theory of strained lattices has been successful 
on the self-organization of elastic stress domains and nanostructures 9l-ll4j|. The elastic 



strain models developed in those studies consider the strain due to a mismatch between two 
different lattices with a coherent interface. Recently, extensive studies have also been made 



on multiferroic nanostructures in epitaxial composite films using phase field model |l5Hl7l|. in 
which a few phenomenological model parameters need to be properly selected for modeling 
the strain field in multiple mismatched lattices. This paper presents a micromechanical 
model based on the theory of elasticity with experimentally measured elastic constants 
and without any fitting paramenter for the strain field in epitaxial nano-composite films 
with impurity nanostructures. In this model, the boundary conditions of the equilibrium 
equation at interfaces are more properly specified as compared with the previous models and, 
therefore, no phenomenological fitting parameter is needed. This model therefore represents 
a considerable improvement from the previous models. While the doped YBCO epitaxial 
film is studied as an example system, the method of the strain analysis developed here can 
be readily generalized to study the self-organization of nanostructures in other epitaxial 
nano-composite films. This paper is organized as follows. In Sec. II, the basic formalism 
of the model is discussed. Studies of the orientation of impurity NRs in epitaxial films on 
lattice-matched and -mismatched substrates are presented in Sec. Ill and IV, respectively. 
Section V contains a conclusion. 

II. Elasticity Model of Impurity Nanostructures in Epitaxial Films 

In this model, the formation of aligned impurity NRs in epitaxial films is assumed to be 
the consequence of the relaxation to the energetically-preferred elastic equilibrium of coher- 
ently strained lattices due to lattice mismatches among film, dopant, and substrate. The 
assumption of the coherently strained lattices is based on the following considerations, (a) 
With a small volume density of dopant, the film and dopant lattices are likely to be coher- 
ently strained as long as the lattice mismatch between film and dopant is sufficiently small. 
Experimentally, semi-coherent interface between YBCO and dopant with some dislocations 



has been observed j4|. The effect of the dislocations is neglected here in order to obtain an 
analytical solution, (b) In the case of the film on a lattice-mismatched substrate, as long 
as the lattice mismatch is sufficiently small, the film could be coherently strained near the 
film-substrate interface until a critical thickness is reached. Beyond this coherent region, 
the strain in the film can be released by introducing dislocations. Since the crystalline con- 
figuration of impurity nanostructures in epitaxial films is strongly influenced by their initial 



formation near the film/substrate interface, only the strain field in the coherent region of 
the film is considered in this study and the effect of non-coherent over-layers is neglected. 
Therefore, the possibility to fabricate impurity NRs with the vertical or horizontal alignment 
in c-oriented films depends on the elastic energy of the strained lattices with respect to other 
possible NR configurations. Note that the horizontal and vertical alignments are the most 
probable, if not only, configurations assuming dopant-film epitaxy is maintained. On the 
(OOl)-cut SrTiOs (STO) substrate, the YBCO ab planes are considered to be twinned. This 
allows simplification of the three-dimensional system to a two-dimensional one that contains 
the [100] and [001] directions. Note that all the calculations can be readily extended to 
the three-dimensional case of orthorhombic lattices. Let Xi and X3 be the components of 
the two-dimensional coordinate along the [100] and [001] direction, respectively. The elastic 
energy of strained lattices can be determined by 

Eei = J EidV+ J E2dV+ J EsdV (1) 

film dopant substrate 

where Ei is the elastic energy density for film (z = 1), dopant (z = 2), and substrate {i = 3), 
respectively. Considering the tetragonal symmetry of twinned YBCO film and the cubic 
symmetry of substrate and dopant, the elastic energy density can be written as jl^ 

Ei = ^Xiiul^ + ^Ai2?i33 + \i3U11U33 + Ai4^i?3 (2) 

where ujk is the strain tensor and Aji = Ci\\ Xi2 = C33, Ajs = 0^3, and Xa = C55 are the elastic 
constants of the material labeled with i. The elastic constants used in this study can be 



found in Refs. 20N23| . Note that the interaction energies at film/substrate and film/dopant 
interfaces are included implicitly in Ei as Ujk is the solution of equilibrium equations with 
the boundary conditions that are the result of the interface interactions. Using the general 
summation rule (also in other formulas in this paper), the equilibrium equations can be 
written as jl^ 

d f dEi \ _^ 
dxk \dujk) 

where ajk = dEi/dujk is the stress tensor. At an interface between two coherently bonded 
lattices, the boundary condition of Eq. ([3]) prescribes continuity of the force on the interface 
and allows for a discontinuity of the strain across the interface. Let n = (rii, n-^) and 
s = (si, S3) be the unit vectors normal and tangential to an interface. The boundary 



condition for the continuity of the force on an interface is 



nfc[(T,fc(l)-a,fc(2)] = 0, (4) 

where o"jfc(l) and (Tjfc(2) are the stress at the interface in lattice 1 and 2, respectively. Along 
the tangential direction of an interface, the lattices on the two sides are deformed in order 
to match each other under the interaction between lattices. As illustrated in Fig. la, if the 
lattice at one side of the interface is stretched (compressed) the lattice at the other side is 
compressed (stretched) accordingly. The total deformation in length of the two lattices along 
the tangential direction is thus assumed to equal the difference of their mismatched lattice 



constants Note that the underlying assumption here is that the bonding between two 
lattices (such as adsorbate-substrate interaction) is much stronger than the lattice bounding 
in each material. This assumption is reasonable for the case of the impurity doped YBCO 
films and many other epitaxial films when lattice mismatch is small or the defects caused 
by lattice mismatch is not a major concern. As the normal components (ukk) of a strain 
tensor describe the change in length per unit length, the discontinuity of the normal strain 
component in the tangential direction of an interface is 

Sk [Mfefc(l) - (1 + fk) Ukk{'2) - fk] = 0, (5) 

where /i = 02/0-1 — 1 and /s = C2/C1 — 1 are the lattice mismatches at the interface along the 
[100] and [001] direction, respectively, and (ai,ci) and (025^2) are the lattice constants at 
each side of the interface. It should be noted that if the above assumption of a strong bonding 
between two mismatched lattices is invalid, the left-hand side of Eq. ([5]) will no longer equal 
to zero. In that case, Eq. ([5]) could be modified by replacing zero with a parameter that has 
to be determined empirically or by model fitting. Moreover, as illustrated in Fig. lb, the 
elastic forces from the two deformed lattices that exert on a surface element perpendicular 
to the tangential direction of the interface should be balanced at the interface, i.e. 

Sk[(Tjk{l)+(Tjk{2)]=0. (6) 

Finally, the strain vanishes deep inside the substrate and the boundary conditions at the 
top surface of a film is simply UkCrjk = 0. It should be noted that in this model the 
lattice mismatches between materials are treated locally at each interface, which allows 
the consideration of different lattice mismatch at different interface. This set of boundary 



conditions is more properly specified than that used in the previous elastic models |l0l.ll5l.ll7]. 
In general, the solution of PDE in Eq. ([3]) is not unique unless its boundary condition is 
properly specified. With the boundary conditions in Eqs. dl])-®, Eq. for epitaxial 
films with impurity NRs becomes a well-posed PDE problem as the solution is unique and, 
therefore, no phenomenological fitting parameter is needed in this model. As far as we 
know this set of the boundary conditions has not been explicitly used before for studying 
the strain field in composite materials. By solving equilibrium equation in each material 
with the boundary conditions, the equilibrium strain for a given configuration of impurity 
nanostructures in a film matrix can be obtained and the elastic energy can then be calculated 
using Eqs. ([T]) and ([2]). In order to test this model, an experimentally well-studied case 



was examined 2J]. In the experiment, a few monolayers of the c-oriented SrLaA104 were 
deposited on STO substrate, where the lattice mismatch between the film and substrate is 
about 4%. The model calculated SrLaA104 lattice deformation of 1.0% expansion in the ab 
)lane and 1.0% contraction along the c axis agrees well with the experimental measurement 



24| of about 1.5% expansion and contraction in the respective directions. 



III. Impurity Nanorods in Epitaxial Films on Lattice-Matched Substrates 

For the case of impurity doping in YBCO films, let's first consider aligned impurity NRs 
in the c-oriented YBCO film on (OOl)-cut STO substrate. To determine the energetically- 
preferred orientation of NRs, the elastic energies for the configurations of NRs aligned in the 
[001] or [100] direction in the film were calculated and compared. Since the lattice mismatch 
between the twinned YBCO ab planes and STO is negligible, only the strain field due to 
the lattice mismatch between film and dopant needs to be considered. In order to solve 
Ujk analytically, the system was further simplified as follows, (a) The length of the NR is 
much larger than its diameter so that the effect of the interfaces between film and dopant 
at the two ends of a NR is negligible. Hence, the strain field is approximately uniform in 
the direction along NRs. (b) The volume density of NRs is small so that the strain decays 
to zero at the halfway point between two neighboring NRs. (c) NRs are assumed to be 
uniformly distributed in the film. Figure 2 illustrates this geometric configuration of NRs 
in epitaxial films for the solution of Eq. ([3]). With these simplifications, Uij was found to 
be linearly dependent on the coordinates for both the configurations of the NR orientation. 
When NRs align in the [001] direction, for example, the equilibrium strain in YBCO was 
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obtained as 



Mil = -AAi3(l -xi/D) 

U33 = - (A11/A13) Mil 
Ul3 = 



(7) 



where xi G [0, D] with xi = at the center of a NR and xi = D a.t the halfway point 
between two neighboring NRs (see Fig. 2a), 



A 



p) [A11W2 + 
AiiAi2 



/3)Al2Wi] 



Wi 



\2 



A,:: 



i2 



(8) 



(9) 



p is the volume density of NRs, and /s and /i are the lattice mismatches between the film 
and dopant in the [001] and [100] direction, respectively. For the twined YBCO a6-planes, 
the a-axis lattice constant in fi is the average of the lattice constants of the a and b axis. 
A similar solution was also obtained for the equilibrium strain in the configuration of NRs 
aligned in the [100] direction. With the obtained equilibrium strains, the elastic energy 
difference between the two orientation configurations was calculated as 



E,i[100] - E,i[001] 



V A?i wi wl /I 



6 [\11W2 + Ai2iyi(l + /s)]' 



1 , ^1 ^ wiAi2(Aii - A12) 
I- P+—p]G p 

W2 J AfiW2 



(10) 

where £'e/[100] and £'e«[001] are the elastic energies of the film with NRs aligned in the [100] 
and [001] directions, respectively, V is the volume of the film, and for /i << 1 and /s << 1, 

G= - V^V-^- (11) 



A12 



Aii(wi + W2) 

Considering p << 1, G as a function of W2 and I/1//3I can be conveniently used as a state 
function for the NR orientation, where W2 is of the elastic constants of dopant. Note that 
< iVi < Aji for a positive-definite elastic energy. When G {w2, I/1//3I) < 0, i^^e/[100] — 
-EejOOl] < and it is not possible to have NRs aligned in the [001] direction and vice 
versa. Hence, G {w2, I/1//3I) = yields a phase boundary that separates the regions in the 
parameter space of {w2, I/1//3I) for dopant, where the vertical or horizontal alignment of 
NRs is not possible in the c-oriented film on a lattice-matched substrate. Figure 3 plots this 
phase diagram for the doped YBCO film where the vertical alignment of NRs is only possible 
in region I. For the BZO or BSO NRs, as shown in the figure, the vertical alignment is the 



energetically preferred state in the c-oriented YBCO 



) nm 



m on a lattice-matched substrate, 
25|. The preference of the NR 



which is consistent with experimental observations 
orientation is determined by the difference in the lattice mismatches between the film and 
dopant in the [001] and [100] directions and the anisotropy of the elastic constants. Because 
/s < /i) the alignment of NRs in the [001] direction results in less deformation in the YBCO 
lattice. Since A12 < An, moreover, the film lattice along the [001] direction is relatively easier 
to deform and, therefore, easier to accommodate NRs. The strain energy is therefore lower 
when NRs aligns in the [001] direction. Recently, efforts have also been made to fabricate the 
YBCO films with the vertically-aligned Y2O3 or CeOs NRs but have not been successful. As 
shown in Fig. 3, Y2O3 and CeOa are all in region II of the phase diagram, which confirms the 
negative result of the experiments. It should be emphasized that the condition obtained here 
is only applicable to the case of the impurity NRs in the film without significant coexistence 
of other forms of impurity inclusions. The existence of secondary impurity inclusions could 
substantially change the strain field in the film and could in turn alter the elastic energy 
minima. This may explain the observed BZO NRs splay with addition of Y2O3 nanoparticles 
in YBCO films [gI,!?]. 

IV. Impurity Nanorods in Epitaxial Films on Lattice-Mismatched Substrates 

In the case of impurity-doped YBCO films on lattice-mismatched substrates, the lat- 
tice mismatch between the film and substrate results in a strained film lattice, which in 
turn changes the lattice mismatch between the film and dopant. Considering that the 
film/substrate mismatch is the same in both the directions of the interface or only in one 
direction, the strain in the film due to the film/substrate lattice mismatch can be easily 
solved from Eq. (|3]) as 

Mil = r(l - xs/h) 

M33 = -^(Ai3/Ai2)mii (12) 

Ui2 = 

where 2:3 G [0, h] with X3 = at the substrate surface, h is the film thickness, 

«;3 + ^'i(l + /.)' ^ ^ 

Wi is given in Eq. 1^, fg = a^/ai — 1, 03 is the unstrained lattice constant of the substrate, 
and ^ = 1 or 2 for the cases of the lattice mismatch in one or two directions of the interface, 
respectively. If fs > 0, i.e. the substrate lattice is bigger than the film lattice, uu > and 
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M33 < 0, and vise versa. The tensile (compressive) strain in the film ab planes due to a 
mismatched substrate leads to a compressive (tensile) strain along the film c axis. 

For the inclusion of impurity NRs in a film on a mismatched substrate, the effect of 
the substrate can be studied approximately by considering the NRs in a pre-strained film 
matrix due to the mismatched substrate. Near the film-substrate interface, the changes of 
the film lattice constants due to the mismatched substrate are 6ai/ai ~ mii(0) = F and 
Sci/ci ~ ^^33(0) = — ^(Ai3/Ai2)F. Consequently, the lattice mismatch between the film and 
dopant becomes /i + 6fi and /3 + 5/3, where 

■5/3 = f (1 + /3)(Al3/Al!)r 

and the state function G{w2, I/1//3I) in Eq. (ITT]) is modified as G{w2, |/i + 5/1I/I/3 + ^/sD- 
Figure 2 plots ^(^2, l/i + ^/il/l/a + ^/sl) as a function of the film/substrate lattice mismatch 
for cases of BZO and BSO in YBCO on an example cubic-lattice substrate that has similar 
elastic constants to STO. It shows that the nanorods align with the c axis (or in the ab 
plane) of the film if fs is smaller (or larger) than a threshold. When fs is larger enough, 
the lattice mismatch between the film and dopant is substantially altered by the strain in 
the film due to the mismatched substrate and, consequently, the energetically-preferred NR 
alignment changes from the [001] to [100] direction. For both cases of BZO and BSO, as 
shown in Fig. 2, this threshold for the transition of the NR orientation is fsc — +1.2% or 
+2.5% for the mismatch in the both or only one lattice direction of the substrate surface. 

The deformation of the film lattice due to the inclusion of the NRs can be calculated 
by averaging the principal components of the equilibrium strain over the film, which can 
be compared with experimental measurement. Because of the different NR orientation, the 
lattice deformation is different in the regions of fs < fsc and fg > fsc- When fg < fg^ the 
deformation of the film lattice calculated with p << 1 is 

Sai ^ Ai3W2(/3 + Sfs) 1 ^ 

ai 2[\iiW2 + il + f3 + Sfs)\i2Wi] 2 

(15) 

fcl ^ AiiW2(/3 + Sfs) Ai3_p 

Cl ~ 2 [XnW2 + (1 + /3 + 6fi)Xi2Wi] 2Ai2 
For fs > 0, 5ai/ai < and Sci/ci > 0, which represents a compression and expansion 
of the film lattice along the a and c axis, respectively. From Eq. (fT5|) . when fs = 0, 
{6ai/ai)/{6ci/ci) = — A13/A11, i.e. the ratio of the lattice-constant changes depends only 



on the elastic constants of the film and is independent of the properties of dopant. In the 
case of BZO NRs in YBCO films on the (OOl)-cut STO substrate, it was measured that 

n 

the YBCO c axis expands about ~ 0.1% ^] while the calculation with Eq. ( 1T5]) yields 
an expansion of 0.2%. The small discrepancy here may be attributed to the presence of 
dislocations around the BZO NRs, which could release the local strain partially 18j. When 
fs > fsc, the deformation of the film lattice becomes 



6ai 1 



(16) 



_W2 + (1 + /i + 5h)wi 

and 5ci/ci = — (Ai3/Ai2)(5ai/ai). In this case, the film lattice is expanded in the ab planes 
and compressed along the c axis. The ratio of the lattice deformations {6ai/ai)/{6ci/ci) 
also depends only on the elastic constants of the film. 

V. Summary and Discussion 

A micromechanical model based on the theory of elasticity has been developed to study 
the configurations of the self-organized impurity nanostructures in epitaxial films. By treat- 
ing lattice mismatch locally at interfaces, the strain field due to multiple mismatched lattices 
of film, dopants, and substrate can be simultaneously considered. Including the effect of 
multiple lattice mismatches is important to the understanding of impurity nanostructures 
in nano-composite films on lattice-mismatched substrates. With a simplified geometry of 
impurity nanorods in YBCO epitaxial films, the strain field in the film was solved in closed 
form. Based on the analytically calculated elastic energy, quantitative criteria have been 
developed in selecting dopants in the YBCO film for possible orientation of the nanorods. 
An excellent agreement between the theoretical predictions and the experimentally observed 
nanostructures has been achieved. The significance of the agreement between the model and 
experiment is that there is no any parameter in the model that can be adjusted to fit the 
calculation to the experiment. The success of this model suggests that the strain field in 
coherent interface layers of lattice-mismatched dopant, film and substrate is the dominant 
driving force for the self-organization of the nanostructures and fine tuning of the strain 
field through engineering selected interface can lead to controllable growth of the desired 
nanostructures in nano-composite films. This micromechanical model can be used as the 
foundation for a numerical study of the impurity nanostructures in presence of dislocations 
and secondary impurity doping in nano-composite films. 

One interesting result of this study is a simple scaling factor between the strains in the 
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directions parallel with and perpendicular to an interface of two mismatched lattices. As it 
can be seen in Eqs. ([7]) and f|T2|) . us^/uu in a strained film lattice depends only on the elastic 
constants of the film and is independent of the lattice mismatch or the elastic properties of 
the other mismatched lattice. In the case of strained epitaxial films on lattice-mismatched 



substrates, this phenomenon has been observed in an experiment [2J]. Experimental con- 
firmation of such a scaling behavior of the lattice deformation in the nano-composite films 
is important and will provide further insights in the strain mediated self-origanization of 
nanostructures in epitaxial films. Such a study could lead to a better understanding and a 
more accurate modeling of the stain field in the nano-composite films. 
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Figure Captions 




Figure 1. Illustration of the boundary conditions of Eq. ([3]) in (a) Eq. ([5]) and (b) Eqs. 
(jl]) and (EI). The open and solid circles are the undeformed and deformed lattices, 
respectively, ai and 02 are the natural lattice constants of the two lattices, Ukk are 
the principal strain along the interface (dashed line), and F„ and Fg are the elastic 
forces along the normal {n ) and tangential (s ) direction of the interface, respectively. 
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Figure 2. Geometric configurations of NRs aligned in (a) the [001] direction (xa-axis) 
and (b) the [100] direction (xi-axis) for solving Eq. ([3]). i? is the radius of the NR 
and 2D is the distance between two neighboring NRs. 
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Figure 3. Threshold of I/1//3I as a function of W2 for the NR orientation in the c- 
oriented YBCO film on lattice-matched substrates. Below (above) the solid curve, the 
vertical (horizontal) alignment is unattainable. The square, triangle, and diamond 
are of BZO, BSO, and Y2O3, respectively. The point for Ce02 is below the bottom 
of the figure. 
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Figure 4. G v.s. the film/substrate lattice mismatch for BZO (a and c) and BSO (b 
and d) NRs in YBCO film on a cubic-lattice substrate with STO elastic constants. 
(a,b) The mismatch is the same in the [100] and [010] directions and the film ab 
planes are twinned. (c,d) The mismatch is only in the [100] direction and the film is 
de twinned. 
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